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ABSTRACT 



Context. Detections of molecular lines, mainly from IU and CO, reveal molecular material in planetary nebulae. Observations of 
a variety of molecules suggest that the molecular composition in these objects differs from that found in interstellar clouds or in 
circumstellar envelopes. The success of the models, which are mostly devoted to explain molecular densities in specific planetary 
nebulae, is still partial, however. 

Aims. The present study aims at identifying the influence of stellar and nebular properties on the molecular composition of planetary 
nebulae by means of chemical models. A comparison of theoretical results with those derived from the observations may provide 
clues to the conditions that favor the presence of a particular molecule. 

Methods. A self-consistent photoionization numerical code was adapted to simulate cold molecular regions beyond the ionized zone. 
The code was used to obtain a grid of models and the resulting column densities are compared with those inferred from observations. 
Results. Our models show that the inclusion of an incident flux of X-rays is required to explain the molecular composition derived 
for planetary nebulae. We also obtain a more accurate relation for the N(CO)/N(H 2 ) ratio in these objects. Molecular masses obtained 
by previous works in the literature were then recalculated, showing that these masses can be underestimated by up to three orders of 
magnitude. We conclude that the problem of the missing mass in planetary nebulae can be solved by a more accurate calculation of 
the molecular mass. 
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1. Introduction 



Molecular material is detected in many planetary nebulae (PNe), 
mainly by the identification of CO and H? lines. CO h as been 
detected in about 100 PNe (Hu ggins et ail 1 19961 |2005l) . while 
infra red emission of H2 has been identified in m ore than 70 
PNe dHora et al]ll999r. ISterling & Dinersteinll2008l) . The obser- 
vations, both in CO and H2, indicate that these objects have in 
general a high or moderate N/O abundance ratio (> 0.3) and 
bipolar morphology. The detections include young and evolved 
PNe, such as NGC 7027 and NGC 7293 (Helix Nebula), sug- 
gesting that molecules can survive over a significant fraction of 
the lifetime of a PN. 

The detection of molecules other than H2 and CO is still re- 
stricted to a few PNe. Most molecules have been detected in 
the well-studied source N GC 7027 (e. g. lLiuet al.|[l996l \l99% 
Hase gawa & Kwokll200lb and in NGC 7293 dTenenbaum et alJ 
2009, and r eferences therein). Fo r other PNe, a leading work 
was done bv lBachiller et al.l (1 1997). They reported line observa- 
tions of some molecules (CO, 13 CO, CN, HCN, HNC, HCO + , 
CS, SiO, SiC2, and HC3N) in a sample of seven objects at dif- 
ferent stages of evolution. Bachil ler et al. I (fl997l) showed that 
the molecular composition in PNe is different from that in stars 
in the asymptotic giant branch (AGB) and from that in proto- 
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planetary nebulae (proto-PNe). As the star evolves from AGB to 
proto-PN and then to the PN phase, CS, SiO, SiC 2 , and HC 3 N 
densities decrease, such that these molecules are not detected 
in PNe. On the other hand, CN, HCN, HNC, and HCO + densi- 
ties increase during the evolution. The HCO + emission is not 
detected or is extremely weak in progenitor objects such as 
IRC+10216 and CR L 2688, but it is relativ ely strong in the PNe 
phase. The stud y of [ Bachiller et alJ (Il997l) i s complemented by 
observations of iJosselin & Bachiller! (|2003), who reported ob- 
servations of 13 CO, CN, HCN HNC and HCO + in a sample of 
compact PNe, and bv lBell et al J (12007). who carried out a search 
for CO + emission in eight objects, including proto-PNe and PNe. 

Black) dl978l) developed the first chemical model for PNe. He 
modeled the ionized gas and predicted the existence of diatomic 
molecules, such as H2, Hj, HeH + , OH, and CH + , in the transi- 
tion zone H + /H°. A steady-state chemical model fo r the neutral 
gas co ntained in globules of PNe was developed bv lHowe et al.l 
(1994). They predicted detectable quantities of CN, HCN, and 
HNC in a carbon-rich gas and concluded that their results agree 
qualitatively with the observations. For HCO + , however, the re- 
sulting density is much lower than that obtained from the obser- 
vational data. 

Based on observations of the ea rly-type PN NG C 7027 by 
the Infrared Space Telescope (ISO). lYan et al.l (IT999) presented 
a thermal-chemical model for the neutral envelope of this ob- 
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ject in a semi-infinite slab approximation. Their work indicates 
that the molecular densities obtai ned for NGC 7027 can be ex- 
plained by a hot gas chemistry. lHasegawa et al.l d2000t) devel- 
oped a spherical symmetric, steady-state chemical model for the 
same PN, with emphasis on radiative transfer, but without a de- 
tailed thermal or dynamical treatment. Their results show that 
the molecular densities in NGC 7027 are a result of a combina- 
tion of photochemistry and high gas temperature. 

iNatta & Hollenbachl (Il998l) developed time-dependent evo- 
lutionary models for the neutral envelope of PNe. Their mod- 
els include the H2 chemistry and thermal balance. They studied 
the effect of shocks, far-ultraviolet radiation, and also soft X-ray 
emission from hot central stars. 

lAli et al.l (12001 1) constructed a time-dependent chemical 
model to investigate the chemistry of clumpy neutral envelopes 
of three PNe: NGC 6781, M4-9, and NGC 7293. Compared 
with data inferred from the observations, the calculated frac- 
tional density of CN is too high by a factor of 2-3, while HCO + 
is less abundant by a factor of 5. Their models also predict too 
much CS and SiO. 

Both INatta & Hollenbachl dl998l) and lAli et al.l d200ll) in- 
cluded X-ray emission from hot stars in their simulations. They 
showed that X-ray photons can aff ect the physical and chemi- 
cal conditions in PNe. In particular. lAli et ail (1200 ll) reproduced 
the density ratios between CO, CN, and HCN by including X- 
ray effects in their simulations. However, observations show that 
the stellar emission is not the only possible source of X-rays 
in PNe . It is now we ll known that X-rays can have a diffuse 
origin dKastnerl 120071 and references therein). One possibility 
is that this emission is produced by a hot bubble, a rarefied 
gas of high temperatur e generated by wind interactions (e.g. 
Mellema & Frank 1995). X-ray emitt ing gas can also be in the 
form of collimated fast winds or jets (Soker & Kastnerll2003l) . 

In this paper we present the results of a self-consistent calcu- 
lation of molecular concentrations in PNe, from the ionized re- 
gion to the external neutral and cold gas. Our calculations span 
a wide range of physical parameters that characterize PNe. In 
order to compare our results with the observational data, we re- 
strict our analysis to the molecules CO, HCO + , CN, HCN, and 
HNC. 

The models are described in Sect. [2] Model results are pre- 
sented and discussed in Sect. [3] while in Sect. [4] our results are 
compared with those derived from observations. A discussion 
about the determination of molecular masses as the missing mass 
problem is presented in Sect. Conclusions are summarized in 
Sect.|6] 



2. Models 

In this section we present the numerical code used to obtain the 
models and the included chemical network. The assumed input 
parameters and the range of values adopted for the free parame- 
ters are discussed. 



2.1. The numerical code 



The numerical code Aangaba dAl eman &Gruenwaldll2011l and 
references therein) simulates the physical conditions in a neb- 
ula illuminated by an ionizing radiation source. The simulation 
starts at the inner border of the nebula and proceeds into the out- 
ward direction. In each position, the physical conditions (e.g., 
atomic, molecular, and electronic densities, and gas tempera- 
ture), as well as continuum and line emissivities, are calculated. 



The thermal and chemical structures are mutually dependent, 
and depend on the incident radiation in the particular position 
in the nebula. Thus, the code performs iterative calculations to 
obtain the solution of the coupled equations. 

The outward-only approximation is adopted for the transfer 
of the primary and diffuse radiation fields. Geometric dilution 
and extinction of the radiation by gas and dust are taken into 
account. 

The density of the gas phase species (atoms, molecules, their 
respective ions, and electrons) are obtained at each position of 
the nebula under the chemical and ionization equilibrium hy- 
potheses. Twelve elements and their ions are included: H, He, C, 
N, O, Mg, Ne, Si, S, Ar, CI, and Fe. 

The original code was improved for also dealing with 
neutral and cold regions. For this, a more extended set of 
molecules and chemical reactions was introduced (details 
in the following section) as well as cooling and heating 
mechanisms expected for a neutral gas. All temperature- 
dependent data, new or previously included in the code, 
were adapted to correctly treat a low-temperature gas. 
The references for data re l ated to the gas tempera t ure ar e 
the following: IXggarwall (fl983l). Aggarwa l et all dl984l): 
Agga rwal! (I1985I) . lHaves & N ussbaumer (1984), Keena n et al 
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(1986), Martin (1988) . ISmitsI d!991h . 
Jl993[> . ISawev & Berringtonl (Il993b IZhang & Pradhan 
dl995l). IStorev & Hummer! (119951). buinet et al l ([199 
Ramsbotto m et al] (1 19971). [Ag garwal & Keenan] dl99' 
CleggetalJ (1 19991) IVemer et al.1 d 19991). iBrav et al.l d 200 
Gupta & Msezand d2000l) . iTavall d2000l) . IWilson & Belli GOO 
and lFerland et al.l d2009l and references therein). 

The gas temperature is calculated under the thermal equi- 
librium hypothesis, i. e. the total gain of energy per unit time 
and volume is balanced by the total loss of energy per unit 
time and volume. Several gas heating and cooling mecha- 
nisms due to atomic species, dust, molecules (H2, CO, OH, 
and H2O), as well as those driven by radiation, are included. 
C osmic-ray heating i s also included assuming the formalism 
oflYusef-ZadehetdJ d2007l) with a n ionization rate of H2 by 
cosmic-ray equal to 1.3 x 10" 17 s" 1 dWoodall et alj|2007l) . Dust 
is included as in Gruenwald et al. (in preparation). Heating and 
cooling mechanisms related t o the H2 molecule are described in 
lAleman & Gruenwald! d20 11 1) . 

The level population of some molecules, required for calcu- 
lating line intensities and the resulting gas cooling, is obtained 
assuming statistical equilibrium for H2, CO, H2O, and OH; local 
thermodynamic equilibrium is assumed for CN, HCN, HNC and 
HCO + . For CO, H2O, and OH the statistical equilibrium were 
carried out using the MOLPOP program developed by M. Elitzur 
(private comm.), which was coupled to the Aangaba code. The 
detailed treatment for the statistical e quilibrium of H2 can be 
found in Alem an & Gruenwald! d201 ll) . The line emissivities of 
CO, H2O, and OH wer e calculated assuming t he esc ape proba- 
bility equation given bv lHollenbach & McKed d!979h . The cool- 
ing due to these molecules is assumed to be the sum of the en- 
ergy emitted by all allowed transitions between the energy levels 
included in the code. Data for transition probabilities, collision 
coefficients, and ene rgy levels were obtai ned from the LAMDA 
dSchoier et al.ll2005l) and HITRAN 2005 dRothman et al.ll2005h 
databases. 



2.2. The chemical network 

Our chemical model includes 95 molecules (see Table [1} be- 
sides the above mentioned atoms and their ions. The chem- 
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ical network consists of 1693 reactions and is based on the 
UDFA 2006 catalog dWoodall et al.ll2007h . The adopted set of 
reactions is appropriate to describe the chemistry of the main 
speci es studied here according to the astrochemical literature 
(e.g. Ivan Dishoeck & Blacklll989t ISternberg & Dalgarnoll 19951: 
Boger & Sternberg 120051) . The rate coefficients are also taken 
from UDFA 2006 (including cosmic-ray ionization rates), except 
those related to photo-processes or to the H2 network reactions. 

The rate coefficients of the photo-processes are obtained 
by integrating the cross-section over the energy distribution 
of the incident radiati on. The cross-sectio ns a re those from 
the Hu ebner database dHuebner et al.l 1 19921) and Ivan Dishoeckl 
d!987l and references therein). When a photo-process cross- 
section is not available, the value of the photoreaction coefficient 
(T) was estimated by a similarity criterion defined as 



a 

r = — r s , 



(1) 



where a (in s _1 ) is the rate in the unshielded interstellar radia- 
tion field given in UDFA 2006 for each process. The R subscript 
denotes "reference", indicating the value of the respective vari- 
able for the process used as reference. We chose the reference 
according to the value of the y parameter: 

- CH photodissociation for 1 .0 > y < 1 .5 

- OH photodissociation for 1.5 > y < 2.0 

- CO photodissociation for y > 2.0 

The y parameters relate the wavelength dependence of the pho- 
toprocess and are also given in UDFA 2006 for each molecule. 

The coefficients for the chemical H2 network (which con- 
tains the H2, HJ, HJ, and H~ species) and details of the dus t 
properties can be found in lAleman & Gruenwaldl (120041 l20lT1) . 
The H2 self-shielding i s incl uded following the formalism of 
iBlack & van Dishoeckl dl987l) . The CO self-shielding and the 
CO cross -shielding by H2 are estimate d from the tabular values 
given bv lvan Dishoeck & Blackl dl988l) . 



Table 1. Molecular species included in the chemical network. 



2 atoms 


3 atoms 


4 atoms 


5 atoms 


6 atoms 


c 2 


NH+ 


C 2 H 


HCO 


CH 3 


CH 4 


CHJ 




NO 


C 2 H + 


HCO + 


CHJ 


CH 4 + 




CH 


NO + 


C 2 N 


HCS 


C 2 H+ 


H 3 CO + 




CH + 


NS 


C 2 N + 


HCS + 


C 3 H + 


H 3 CS + 




CN 


NS + 


c 3 


HNC 


C 2 NH + 


NH+ 




CN + 


o 2 


c + 


NH 2 


H 2 CO 






CO 


o 2 + 


CH 2 


NHJ 


H 2 CO + 






co + 


OH 


CHJ 


2 H + 


H,CS 






cs 


OH + 


CNC + 


OCS 


H 2 CS + 






cs + 


SiC 


C0 2 


SiH+ 


H 2 NC + 








SiC + 


CO+ 


SiOH + 


H 3 + 






H 2 + 


SiH 


H 2 


so 2 


H 3 S + 






HeH + 


SiH + 


H 2 + 




HCNH+ 






HS 


SiO 


H 2 S 




HCO} 






HS + 


SiO + 


H 2 S + 




NH 3 






N 2 


SO 


H 3 




NH+ 






N 2 + 


SO + 


HCN 




SiCH+ 






NH 




HCN + 











2.3. Parameters of the models 

We obtained a grid of theoretical models with input parameters 
in ranges typical of PNe to the understanding of the correlation 
between molecular concentrations and some properties of PNe. 

The primary energy sources are the central star and the X- 
ray emission produced in a hot bubble located in a cavity be- 
tween the central star and the main PN shell. The central star is 
assumed to emit as a blackbody, defined by its temperature (T.) 
and luminosity (L„). The hot bubble emission is represented by a 
central source of X-rays, defined by the spectral distribution (0.3 
- 2 keV) and by the integrated X-ray luminosity (Lx). The X-ray 
spectral distrib ution assumed is that derived for NGC 5315 by 
lKastnejd2007h . 

The gas density and chemical composition, as well as the 
dust properties (density, composition and size) describe the neb- 
ula. The nebular gas density is represented by the total number 
of hydrogen nuclei («#)■ Most models were obtained with a ho- 
mogeneous density. For testing the effect of density profiles in 
the molecular composition, models with different density pro- 
files are also analyzed. The tested radial profiles are the follow- 
ing: a homogeneous distribution, a power law (r~ p , 1 < p < 4), 
a distribution defined by equilibrium pressure, or a combination 
of them. 

The nebula is assumed to be spherically symmetric. The in- 
ternal radius of the nebula is defined as Rq = 10 15 cm. The exact 
value of Rq does not affect the results presented in this paper for 
values up to 30 % of the ionized radius. 

The adopted chemical composition, homogeneous through- 
out th e nebula, corresponds to the mean values for PNe accord- 
ing to lKingsburgh & Barlowldl994l) . The abundances of Mg, Si, 
CI, and Fe are not provide d by these authors. In t hese c ases, we 
adopted the values given by Stasiri ska & Tvlendald 19861) to make 
a rough correction for grains depletion. Since these elements do 
not have much effect on the gas cooling or the resulting molec- 
ular densities, their exact proportion in the form of grains is not 
important. The abundances of C, N, and O may affect the result- 
ing molecular composition. Therefore, we also analyzed models 
with different values for the abundance of these elements. 

Dust is included as graphite spheres of radius 10~ 2 //m. 
The dust-to-gas mass ratio (Md/M g ) is assumed to be constant 
throughout the nebula. 

The adopted ranges for the free parameters for the planetary 
nebulae models are given in Table|2] with the corresponding ref- 
erences. A standard model with a given set of input parameters 
was chosen, with the corresponding adopted values also listed 
in Table [2] To study the effects of each parameter individually, 
we varied one of the parameters within its typical range, while 
keeping the others fixed. Hereafter, unless otherwise noted, the 
parameters of the models are those of the standard model. 

In the following sections the discussion of the molecular 
composition is made mainly in terms of the column density of 
the CO molecule, N(CO), since this value is easily obtained from 
the observations, where the column density is the integral along 
the radius outward from the central star (both hemispheres are 
taken into account). Models are then obtained with the calcula- 
tions stopping at different values of N(CO). 



3. Model results 

We outline here the thermal and chemical structures for the stan- 
dard model. This first analysis provides a useful guide for the 
discussion about the molecular chemistry in models with differ- 



3 



Kimura, Gruenwald and Aleman: Molecular chemistry and the missing mass problem in PNe 



Table 2. Model parameters: variation range, references and stan- 
dard value. 



Parameter 


Variation range 


References 


Standard value 


r,(K) 


3 x 10 4 - 3 x 10 5 


1, 2 


10 5 1 


U (L G ) 


10 2 - 1.2 x 10 4 


1,3 


3 x 10 3 1 


n H (cm" 3 ) 


10 2 - 5 x 10 5 


4, 5, 6 


10 5 1 


C/H 


10~ 4 - 6 x 1(T 4 


7, 8 


5.50 x 10~ 4 1 


N/H 


io~ 4 - 6 x 1(T 4 


7, 8 


2.24 x 10~ 4 1 


O/H 


io~ 4 - 6 x itr 4 


7, 8 


4.79 x 10~ 4 1 


M d /M g 


io~ 3 - 10~ 2 


9 


5 x 10~ 3 1 


L x (ergs" 1 ) 


0; 10 30 - 10 32 


10, 11 


5 x 10 3 ' t 



...es: (1) Bloecker (1995); (2) Phillips (2003); (3) Phillips (2005); (4) 
lLiuetal.1 JlOOlh; (5) ISt anghellini & KaleJ 4l989l)~(6) iKingsburgh & E nglish 
119921); ml Kingsburgh & Barlow! 419941): (8) iMilanova & Kholtyginl 420091): 
(9) IStasinska & Szczerbal 419991) : < 10) iKastneJ 42007b : ( 1 1) iMellema & Frankl 
419951) . 



ent free parameters and for a comparison with the observations, 
which will both be presented in Sect. [4] 



3.1. Basic chemical structure of PNe 

Since the central star of PNe is an intense source of UV photons, 
ionized matter is present in any PNe. If there is enough mass to 
absorb the ionizing and dissociating photons, a given nebula can 
also have an outer region of neutral and cold matter (radiation- 
or ionization-bounded nebula). 

To simplify our discussion about the molecular chemistry in 
PNe we can identify some regions according to the predomi- 
nant form of H, C, or O (ionic, neutral or molecular). We identi- 
fied the following regions: H + , H°, H2, and CO. As the name 
suggests, the main characteristics of the three first regions is 
the predominant form of hydrogen. The CO region is defined 
as the region where the CO molecule locks up all oxygen or 
carbon (the less abundant), i. e., where CO is fully associated. 
With this definition, the CO region corresponds to the external 
and coldest zones where hydrogen is molecular. We freely de- 
fined boundaries for these regions; from the inner to the more 
external regions, the boundaries are defined by the conditions 
n(H + ) = n(H°), n(H°) = n(H 2 ), and n(CO) = n(X° ), where 
n(X) is the volumetric density of the species X, and Xc,o repre- 
sents the less abundant element (C or O). Indeed, the transition 
from one region to another is smooth and the extent of the tran- 
sition zones depends on the stellar and nebular characteristics. 
The transition regions may be relevant for the molecular produc- 
tion, since the coexistence of two or more forms of a given ele- 
ment may enhance or inhibit some molecular processes. For the 
present discussion, we distinguish three transition zones here: 
H + /H°, H°/H 2 , and C + /C°/CO. 

The physical conditions in the H + region are regulated by 
energetic photons emitted by the hot central star (T t > 3 x 10 4 
K). In this region, electrons are produced mainly by the ioniza- 
tion of hydrogen and helium. The gas temperature is typically 
about 10 4 K. In this harsh environment molecules can survive 
just in the transition zone toward the H° region. In this transition 
zone (H + /H°) the coexistence of cations, anions, and the neutral 
atom of hydrogen in a warm temperature environment provides 
the formation of H 2 by alternative routes and not just by grain 
surface reactions (Alem an & Grue nwald 2004|). 



In the H° region hydrogen is predominantly atomic. Here, 
the UV photons above 13.6 eV are absorbed by gas and dust in 
inner shells, while the flux that dissociates H 2 is intense. The 
physical conditions of the gas in this region are controlled by 
■far-UV photons (6.0 - 13.6 eV) with high influence of X-rays 
photons (0.3 - 2.0 keV). The ejection of electrons from grains 
(photoelectric effect) is the main source of heating, balanced by 
cooling due to metals (mainly C, N, and O). The electrons are 
provided mainly by the ionization of H° (n(H + )/n# ~ 10~ 4 ), 
maintained by X-ray photons, and from the C° ionization by far- 
UV photons. Nitrogen and oxygen are predominantly in atomic 
form. This region extends up to the region where photons able to 
dissociate the molecular hydrogen are mostly absorbed. Under 
this condition, the hydrogen becomes predominantly molecular, 
characterizing the H 2 region. 

In the H 2 region the heating and most chemical reactions are 
dominated by radiation. The cooling by CO becomes significant 
and dominates the total energy loss. Hydrogen is predominantly 
molecular (H 2 ) and has a strong influence on the chemical re- 
actions, since H 2 facilitates the formation or hydrogenation of 
some molecules. Nitrogen and oxygen are still in atomic form. 

Depending on the spectral distribution and intensity of the 
incident spectrum, as well as on the nebular parameters, a neb- 
ula can be extended enough such that the radiation able to ionize 
C° or to dissociate CO is attenuated in the more external region. 
This is what we call the CO region, where the CO dissociation 
rate is very low. In this region the main source of energy grad- 
ually changes for increasing radius from the photoelectric effect 
to heating by cosmic rays. The cooling due to the CO molecule 
is the dominant mechanism of energy loss, and the electron den- 
sity is determined by the ionization of C° or of heavy metals of 
low ionization potential (Mg, Fe), where C + abundance is low. 
Nitrogen is predominantly in molecular form (N 2 ), while oxygen 
is in CO form. In this region the temperature reaches its lowest 
value (T * 40 K). 

3.2. Main molecular formation and destruction routes 

In this section we analyze the molecular distribution along the 
nebular radius. The molecular distribution is discussed in terms 
of the volumetric density (cm 4 ), since this quantity provides lo- 
cal information about the molecular processes (and consequently 
the physical conditions) in different regions of the nebula. 

Figure[TJshows the gas temperature and relative (to n#) abun- 
dances of C°, C + , N°, N 2 , CO, HCO + , CN, HCN, and HNC 
versus the CO column density, N(CO), for the standard model. 
The vertical lines in the figure indicate the boundaries of some 
of the previously defined regions. Following our definition, the 
lines corresponding to the standard model (from left to right) are 
n(H°) = n(H 2 ) and n(CO) = «(O ). The regions showed in the 
figure correspond, from left to right, to the H°, H 2 , and CO re- 
gions. The H + region is not included in Fig. Q] since molecular 
concentrations are negligible in this region. The dependence of 
N(CO) on the radius for the standard model is shown in Fig. [2] 

To understand the pattern presented by the density of the 
molecules along the nebula, shown in Fig. Q] we identify the 
main formation and destruction processes. 

Since the reaction chains of the molecules CO and HCO + 
are intrinsically coupled, their densities are mutually dependent. 
Figure [3] presents a schematic diagram of the main processes 
between CO and HCO + . Details of the reactions summarized in 
the diagram are given below. 
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Fig. 1. Atomic and molecular abundances and gas temperature 
versus N(CO) for the standard model. The line at the top indi- 
cates the H , H2 and CO regions. (A color version of this figure 
is available in the online journal). 
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Fig. 3. Schematic diagram showing the chemical paths for for- 
mation and destruction of CO and HCO + . 



One path for the formation of the radical HCO + is via the 
reaction 



CO -4 HCO + . 

This reaction is efficient only in the CO region because it re- 
quires a high survival rate of H3 , which is expected to be very 
short-lived because of its rapid dissociative recombination in the 
presence of electrons, a consequence of an ultraviolet radiation 
field. An X-ray incident flux raises the density of HJ by increas- 
ing its production rate; nonetheless, this is not the main source of 
HCO + if there is strong UV radiation. In the H° and H2 regions 
where the UV radiation field is intense, the main formation route 
ofHCO + is 

CO + 4 HCO + . 

The molecular ion CO + is formed mainly by the route 



o + ^ OH + 4 



H 2 + ^ H 3 + 



OH ^ CO + . 



Fig. 2. Dependence of N(CO) on radius for the standard model. 



An X-ray incident flux is required for the formation of CO + 
molecules, since it contributes to the increase of the electron 
and + densities. The most favorable condition for the CO + for- 
mation is found in the transition zone H°/H2, where there are 
enough H2 molecules for the synthesis of OH, but their density 
is not high enough to efficiently convert CO + into HCO + . In the 
CO region the main formation route of CO + is the ionization of 
CO by cosmic rays at a low production rate, leading to a very low 
density compared to the CO + formed in the H° and H2 regions. 

CO is formed via HCO + + e -> CO + H and via CO + + 
H° — > CO + H + , although the contribution from the latter is 
minor. The destruction of CO is dominated by radiation, not only 
by dissociation, but also by reactions involving He + , which is a 
product of ionization by X-rays. 

The scenario is more complicated for CN, since there are 
many reactions that can play the major role in forming the 
molecule across the nebula. This can be verified in Fig. |4j which 
presents the formation rate of the main reactions leading to the 
CN molecule as a function of N(CO). 

For HCN, the main formation reaction in most of the nebula 

is 



HCNH + 4HCN. 
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Fig. 4. Rates of the main reactions forming CN for the standard 
model. (1) H° + CN + CN + H + ; (2) HCN + + e~ -> CN + H°; 
(3) HCN + photon -> CN + H°; (4) HCNH+ + e" -> CN + H 2 ; 
(5) C° + NO -> CN + O ; (6) CH + N° -> CN + H°. 



HCNH + can be produced by different routes triggered by ions 
generated by photons or by cosmic rays. The main route de- 
pends on the specific region and can include the reaction with 
activation barrier C + + H2 — » CH + +H°, which is processed in 
the H°/H + transition zone (where T gas > 10 3 K). The charge 
exchange between H° and HCN + also contributes in the H° re- 
gion and C° + NH2 with some importance in the transition zone 
C + /C°/CO. Photodissociation is the main destruction reaction in 
the H° and H2 regions, as in the inner zone of the CO region. 
When photons able to dissociate HCN are scarce, the destruc- 
tion is dominated by collisions with products and subproducts 
of the ionization by cosmic rays. 

The destruction of CN and HCN is dominated by photodis- 
sociation in most of the nebula. In outer zones of the CO region, 
where the photodissociation is not effective, the destruction is 
controlled by collisions with products and subproducts of ion- 
ization by cosmic rays, such as H3 and HCO + . 

The chemistry of HNC can be described just by substituting 
HCN by HNC, with small differences in the rate coefficients of 
the reactions. 



3.3. The effect of C/O and N/O ratios on the chemistry 

To discuss the C/O ratio we can distinguish two regimes for the 
molecular production: one in which the radiation controls the 
chemistry (in the H° and H2 regions) and another in the CO re- 
gion, where the main chemical processes are induced by cosmic- 
rays. In the H° and H 2 regions the molecular chemistry is not 
very sensitive to the C/O ratio, such that the molecular concen- 
tration in an oxygen-rich gas and in a carbon-rich gas is similar. 

On the other hand, CO exhausts the less abundant element 
(C or O) in the CO region, such that the molecular densities 
in a carbon-rich gas are significantly different from those in an 
oxygen-rich gas. C-bearing molecules, such as those of the CH„ 
and C„ families, are highly favored in a carbon-rich gas, while 
OH, H2O, and NO are abundant in an oxygen-rich gas. The 
chemistry of second-row elements, such as Si and S, is also only 
affected by the C/O ratio in the CO region. In the specific case of 
the species CN, HCN, and HNC, in the CO region their presence 



is highly favored in a carbon-rich nebula, which does not occur 
if the gas is oxygen-rich. Figure [5]illustrates these statements. 

We also studied models with different N/O ratios, since PNe 
detec ted in CO and H2 have, in general, hig h or moderate N/O 
("e.g. lHuggins et alJll996tlKastner et alj|1996h . Our results show 
that models with different nitrogen abundances result in similar 
molecular concentrations. Thus, the N/O ratio correlation with 
the detection of molecular lines in PNe may be mainly related to 
other properties of PNe that present moderate or high N/O ratio 
(for example, high progenitor mass, high central star tempera- 
ture) rather than the abundance of N by itself. 

4. The importance of X-rays for the molecular 
chemistry 

In this section we compare our results with observational data 
from the literature. The results are presented and discussed in 
terms of column densities, since molecular data obtained from 
the observations are generally given by these quantities. 

In Fig. [6] we compare column density ratios inferred from 
observations with our results. Ratios are given as a function of 
the column density o f CO. Symbols repres e nt observational data 
from the literature (Bachill er et al.l 119971 Lfosselin & Bachilierl 
I2003h and curves correspond to our models. Two of these mod- 
els include X-ray emission from a hot bubble, with different 
values for Lx- Two models for which the central star is the 
only energy source are also shown: one model with a very hot 
star (T, = 3 x 10 5 K) and another with a low-temperature star 
(T* = 5 x 10 4 K). Other parameters are the same as those from 
thestandard model, except for the dust-to-gas mass ratio which is 
equal 1(T 2 . As can be seen in the figure, the A^(HCO + )/A^(CO) 
and N{CN)/N{CO) ratios are not reproduced by models with- 
out significant quantities of X-rays (situation represented by the 
model with a low-temperature central star). 

The molecular chemistry in PNe is sensitive to the stel- 
lar temperature mainly because hot stars can emit a significant 
quantity of X-ray photons (hv > 50 eV). Owing to the low 
ionization cross-sections for X-rays, these photons can cross 
long distances in the nebula. As a result, they can increase the 
ionization rate in regions where the gas is predominantly neu- 
tral, favoring ion-molecule reactions. Additionally, X-ray pho- 
tons can maintain a relatively high electron density in those re- 
gions, which can determine molecular density ratios. The im- 
portance of the effect of X-rays on the chemistry was stressed by 
studies of XDRs, an acronym for X-ray dominated regions (e.g. 
iMalonev et alll996l) . 

As can be seen in Fig. [6] the model with T, = 3 x 10 5 K 
(and L x = 0) shows that for the N(HCO + )/JV(CO) ratio, the data 
can be explained by models that include a very hot central star, 
source of significant X-ray emission. In the set of models with 
Lx = 0, the minimum temperature of the central star that could 
fit the data for N(UCO + )/N(CO) is T* ~ 1.5 X 10 5 K. 

However, notice that even the lowest values of 
N(CN)/N(CO) cannot be reproduced by models that do 
not include hot bubble X-ray emission. Indeed, if a hot bubble 
emission is taken into account, at least in the range of values 
inferred by the observations (10 31 erg s _1 < Lx < 10 32 erg 
s ), the effects of X-rays emitted by the star are comparatively 
negligible. 

Theoretical values of the JV(HNC)/iV(HCN) ratio only re- 
produce the highest values obtained from the observations. Since 
HCN and HNC have similar chemical paths, with similar rate co- 
efficients, the value of iV(HNC)/iV(HCN) is near unity in most 
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Fig. 5. Relative abundances for a carbon-rich gas (C/O = 1.5, solid lines) and for an oxygen-rich gas (C/O = 0.66, dashed lines) 
for: (a) O-bearing molecules; (b) C-bearing molecules; (c) S-bearing molecules; (d) CN and HCN. The vertical line indicates the 
boundary between the H2 and CO regions. The line at the top indicates the H2 and CO regions in the nebula. 



of our models. The lowest observational values are not repro- 
duced by any set of free parameters. In the astrochemical lit- 
erature, the derived 7V(HNC)//V(HCN) from observations in the 
interstellar medium ranges from 0.01 to 10 dTurner et al.lll997t 
lHirotaetal.ll 19981: IStauber et al.ll2004i) . This range of values is 
only poorly understood yet. Isomeration, alternative formation 
routes, or reactions that are processed efficiently for one of 
those molecules and not for the other in some particular phys- 
ical condition are some possible explanations proposed in the 
literature (eg. iGoldsmifh et all 1 19811 119861: iTurneret al.lfl997t 
iHirota et al.lfl998l;IStauber et al.ll2004l) . 

The high 7V(HCN)//V(CN) ratio compared to the obser- 
vational data indicates that TV(HCN) is underestimated, since 
7V(CN) is well-reproduced compared to N{CO) in our mod- 
els. The difficulty in reproducing the MHCNVMCN) ratio 
was also found bv powe et alJ (U994I) lAli etaLl (l200ll) . and 
iRedman et alj d2003 )~ lHasegawa et al. d2000h only reproduced 
the low value of A^(CN)/AT(HCN) for NGC 7027 by postulating 
a high-temperature gas in the neutral region (T = 800 K), where 

the conversion reaction CN — >HCN is efficiently processed. In 
our models, the temperature distribution is consistently calcu- 
lated assuming thermal equilibrium. Our models do not show 
neutral regions with such high temperature and, therefore, the 



reaction above is not efficient. Another heating mechanism (as 
shock for example) should therefore be present. Another possi- 
bility that could explain the low HCN density in the models is 
the a bsence of an important reaction in the chemical database 
(e.g. iTurner et al.|[l997l; IStauber et al.ll2004l) or inaccurate rate 
coefficients. 

In brief, most of the data inferr ed from t he ob servations 
can be explained by X-ray radiation. lAli et all (12001 ) also con- 
cluded that X-rays dominate the chemistry and lead to the ob- 
served molecula r compos i tion in the three PNe studied by them. 
However, while lAli et al.l (1200 ll) justified the existence of high 
concentrations of HCO + by the presence of H3 , our results show 
that in most cases (N{CO) < 10 17 crrT 2 ) the main formation 
processes occur via CO + , whose density is also enhanced by X- 
ray radiation. Moreover, lAli et al.l (12001 1) justified the inclusion 
of X-ray emission in their models because of the high temper- 
ature of the central stars. Although the X-ray emission from a 
hot star is a possible explanation for 7V(HCO + )/7V(CO) ratio, our 
study shows that the 7V(CN)/7V(CO) ratio is only adequately re- 
produced by models including X-ray emission from a hot bub- 
ble. Moreover, stellar temperatures would need to be atypically 
high to reproduce the 7V(HCO + )/7V(CO) ratio inferred from the 
observations. 
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Fig. 6. Column density ratios versus N(CO). Lines correspond to models T* = 5 x 10 4 K and = erg s _1 (solid lines); T» = 
10 5 K and L x = erg s 1 (dashed lines); T» = 3 x 10 5 K and L x = erg s 1 (dotted lines); T. = 10 5 K and L x = 10 31 erg s 1 
(dot-dashed lines); T* = 10 5 K and L x = 10 32 erg s~ 1 (dot-dot-dashed lines). Other parameters are those from the standard model. 
Squares and triangles represent observational data, with triangles indicating upper or lower limits, according to their orientation. (A 
color version of this figure is available in the online journal). 



As said above, the curves presented in Fig. [6] correspond 
to models with input parameters of the standard model, except 
those for the stellar temperature and for the X-ray flux. Some 
deviation is expected for different sets of parameters. As an ex- 
ample, in Fig.H the N(HCO + )/N{CO) and JV(CN)/W(CO) ratios 
are presented for models with L„, and Mj/Mg in the range 
given in Table [2] and L x and T m the same as for the standard 
model. The parameter for which the molecular chemistry is more 
sensitive, besides the X-rays discussed above, is the dust-to-gas 
mass ratio (the highest and lowest curves in the figure for both 
ratios are for models with M d /M g = 10~ 2 and M d /M g = 1(T 3 , 
respectively). The presence of dust affects the molecular chem- 
istry in three ways: (1) through absorption of the radiation, (2) 
because they act as catalysis for the formation of H2 molecules, 
and (3) through the photoelectric effect, which heats the gas. 
As a consequence, for higher dust-to-gas ratios the H2 region 
is wider, the H° region is narrower, and both regions are closer 
to the radiation central source. Consequently, the density of ion- 
ized species is higher and the ion-molecule reactions are more 
effective. Although the effect of dust on molecular chemistry is 
significant, it is not possible to reproduce most of the column 



density ratios inferred from the observations without including a 
source of X-rays. 



The molecular composition, as well as the physical con- 
ditions in each region (gas temperature, ionic distribution), in 
a plot N(X)/N(Y) versus N(CO) is weakly dependent on the 
gas density and on the stellar luminosity. Accordingly, differ- 
ent density profiles do not change the discussion about the ra- 
tios between molecular densities. The only exception is the ratio 
N(CN)/N{CO). Since the concentration of these two molecules 
are not strongly connected, small differences in gas temperature, 
ionic abundances, and incident flux convert in non-negligible 
differences in the referred ratio, mainly in the H2 region. We em- 
phasize that this ratio is sensitive to n# and L„ although the ef- 
fects of the other parameters (X-ray luminosity, star temperature 
and dust-to-gas mass ratio) discussed above are more significant. 
In general, for higher gas density (or lower stellar luminosity), 
the ratio N(CN) /iV(CO) is higher. 
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Fig. 7. N(HCO + )/N(CO) (top panel) and N(CN)/N(CO) (bot- 
tom panel) versus N(CO). The gray area indicates the range of 
values obtained with our grid of models with L*, n# and Md/M g 
in the range given in Table [2] and Lx and T» the same for the 
standard model. The solid curve represents the standard model. 
Squares represent observational data. 



5. The molecular mass and the missing mass 
problem in PNe 

Masses of PNe progenitors in the main-sequence can be as high 
as 8 M as i ndicated by the detection of white dwarfs in open 
clusters (e.g. lKalirai e"tal]|2007l 120081: 1 Williams et al.ll2009l) . On 
the other hand, masses obtained for the central stars of PNe and 
for the nebular component c an be as high as 1.2 and 0.3 M , re- 
spectively (e.g. Kwok 2000). Consequently, there must be some 
mass that is not being accounted for. The difference between the 
total mass of a PNe and that of the main-sequenc e progenitor i s 
known as the missing mass problem in PNe (e.g. lKwokl [l994). 
Because the above calculation of the n ebular masses only takes 
into account the ionized mass, iKwokl (1994) suggests that the 
missing mass problem could be solved if the neutral matter in 
PNe is taken into account. 

A more recent determina tion of the total mass as s ociate d 
with PNe was obtained by iBernard-Salas & Tielensl (|2005), 



who calculated for two objects the sum of the nebular mass 
(ionized, atomic, and molecular) and the mass of the central 
sta r. Using molecular masses derived from CO observ ations 
bv lHuggins etal] dl996l) . IBernard-Salas & Tielensl (120051) found 
good agreement between the mass expected for th e progenitor 
star of NGC 6302 (~ 4.5 M , iMarigo et al.ll2003h and the to- 
tal mass estimated for this PN (~ 3.9 M ). On the other hand, 
for NGC 7027 part of the ejected mass seems to be missing 
(~ 2 - 3 M ), even taking the neutral matter into account. It must 
be noted that for the former PN the atomic mass corresponds to 
almost the whole nebular mass, while for the second one most 
of the gas is molecular. 

The molecular masses derived by Huggin s et al.l (1 19961) are 
obtained from the intensity of the rotational line 2-1 of the CO 
molecule and approximations for the ratio between the CO and 
H2 column densities. Fro m our models, howeve r, we conclude 
that the method adopted bv lHuggins et alJ (|1996) can underesti- 
mate the molecular content of PNe, as is di scussed below . 

The molecular mass was obtained by Huggins et al. (1996) 
from the rotational line 2-1 of CO molecule according to the 
expression 



M, 



mol 



2.6 X lQ- w FD 2 /f, 



(2) 



where F is the CO line flux (in K km s _1 arcsec 2 ), D is the dis- 
tance to the PN (in kpc), / is the column density of CO rela- 
tive to that of hydrogen, and M mo \ is given in solar masses. The 
formula includes a correction for the presence of helium with a 

representative abundance val ue of He/H — .1. 

To obtain a value for /, iHuggins et al.l (1 19961) assumed the 
following approximations: (1) all hydrogen nuclei are in the 
molecular form; (2) CO is fully associated. With these simplifi- 
cations the authors adopted a value for / equal to 3 x 10 4 , based 
on the C/H and O/H ratios in PNe where both abundances have 
been measured. Then, for / 



/ 



jV(CO) _ jV(CO) 
JV(H) ~ 2N(H 2 ) 



3 x 10 



-4 



(3) 



Our models, however, show that CO is fully associated only 
in the CO region (see Fig.[TJ, which is mostly dissociated in the 
H2 region. As a consequence, the ratio between integrated quan- 
tities, such as column densities, depends on the specific nebula, i. 
e., on its total mass and how the molecules are distributed along 
the nebula. As stated above, the distribution of the molecules 
depends, in turn, on the incident radiation spectrum and on the 
nebula parameters (density, elemental abundance, etc.). 

Figure shows the N(U 2 )/N (CO) ratio (= 0.5/" 1 ) for a set 
of models that spans the range of physical parameters given in 
Table|2](excepting models with a chemical composition different 
from that adopted for the standard model and models with Lx = 
0). 

Based on our models we can represent the N(H2)/N(CO) ra- 
tio according to the equation 



log ^£ = -107.5 +22.507 C o- 1.4537^ 



JV(CO) 



CO 



+3.028 x 10~ 2 y£ o , 



(4) 



logJV(CO). 

-2 



This relation is valid for 10 < 



where Yqo 

N(CO) < 10 19 crrT 2 and is derived for the standard model. It 
corresponds to the solid curve in Fig. [8] For the other models 
shown in this figure the maximum deviation is of a factor two. 
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Fig. 8. Column density ratio N(H2)/N(CO). The gray area indi- 
cates the range of values obtained with our grid of models. The 
solid line represents the values obtained for the standard model 
(Eq.|U). T he horizon t al das hed line corresponds to the hypothe- 
sis o f THuggins et all d!996l) that CO is fully associated, applied 
to the chemical composition of the standard model. 



For models with Lx = 0, not shown in the figure, the results can 
differ by up to a factor five. 

Taking into account the results from models with values for 
the C and/or O abundance different from that of the standard 
model, the N{H2)/N(CO) ratio derived from Eq. |4]must be mul- 
tiplied by the following correction factor 

4.79 x 1(T 4 



Xqo/H 



where XqoIH is the abundance of the less abundant element 
between C and O. 

As said above, IHuggins et alJ (1 19961) assumed that CO is 
fully associated. The horizontal dashed line in Fig.|8]corresponds 
to this hypothesis applied to the chemical composition of the 
standard model. Note that the assumption is only correct for very 
high values of the CO column density, as can be seen in Fig. [8] 

The distribution of PNe with CO column densities calculated 
by Huggins et al. (1996) is shown in Fig. [9] Comparing this his- 
togram with Fig. [8] it can be noticed that there is a large number 
of PNe in their sample for which their approximation for the 
ratio N(H2)/N(CO) is not valid. For these objects the obtained 
molecular mass is then underestimated. 

It is important to keep in mind that when IHuggins et al.l 
( 1996) determined the column density, they did not take into ac- 
count individual characteristics of each nebula, such as geometry 
and angular size. That is, they assumed a simplified correction 
factor to take into account the beam dilution for objects with di- 
ameters smaller than the beam. In these cases, the observational 
CO column density can differ from the radial column density, as 
obtained from the models. 

The molecular masses obtained by Huggins et al. (1996) 
for their sample can then be recalculated from Eq. |2] assum- 
ing the JV(H 2 )/W(CO) ratio (= 0.5/" 1 ) from Eq. g] scaled to 
Xqo/H = 3 x 1CT 4 implicitly assumed b y them. The new val- 
ues for the masses and those obtained by Hugg ins et al.l (fl996) 
are compared in Fig. [TO] As expected, the differences between 



log N(CO) (cm 2 ) 



Fig. 9. Distribution of N(CO) for PNe. Data from Hugginsetaj] 
(1996). 



new and previous masses are greater for nebulae with lower val- 
ues for N(CO). Our results do not show the tren d of increasing 
molecular mass with increasing N(CO) noticed by Hu ggins et al.l 
d!996l) . It must be remarked that N(CO) is not necessarily related 
to the nebular size and mass, but depends on the distribution of 
the CO molecule on the nebula which, as said above, depends 
on the incident radiation and on the absorption properties of the 
nebula. 
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Fig. 10. Molecular mass vers us the CO column den sity. Squares 
show the masses obtained by iHuggins et alJ dl996l) . and circles 
indicate masses obtained by the method described in the text. 



According to observations and evolutionary models, masses 
for the central star of PNe are mo stly in the range 0.55 M to 
0.65 M Q (e.g. Stas iriska et alJ[l997h . which corresponds to pro- 
genitor stars with m asses approximately in the range 1 M to 3 
M Q (e.g. lKwokll2000l) . Thus, the total nebular mass (given by the 
difference between the mass of the progenitor star and that of the 
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central star of the planetary nebulae) should be in the range 0.35 
M to 2.45 M . 

From masses for the ionized gas derived fro m radio obser- 
vation s and their calculated molecular masses, Huggin s"et al.l 
(1996) obtained values for nebular masses in the range 10 3 to 
1 M . These results do not reproduce the range given above for 
the expected nebul ar masses of PNe. Th e distribution of nebu- 
lar masses given by Hug gins et al.l ( 1 19961) is shown in the upper 
panel of Fig.Qj] 

The nebular mass can also be calculated from our mod- 
els. However, to obt ain the total mass for the nebulae of the 
Hugg ins et alj (1 19961) sample, a specific model for each nebula 
should be obtained. This is beyond the scope of this paper since 
we merely aimed for an estimate of the nebular mass, after the 
correction of molecular masses propo sed in the present p aper. 
Summing the ionized mass adopted bv lHuggins et alJ d 19961) and 
the corresponding molecular mass recalculated as shown above 
for each planetary of their sample, the obtained nebular mass 
distribution is that shown in the lower panel of Fig. Q~T] The new 
calculated masses are up to 4.5 M Q , which agrees better with the 
expected distribution of nebular masses. Note that this is still a 
lower limit, since the atomic mass is not taken into account. 




Fig. 11. Nebular mass distribution for the dHuggins et al.l fl996) 
sample: their calculation (top panel), and after recalculation of 
the molecular mass according to the method described in the 
text (botton). 



6. Summary 

A self-consistent numerical code was adapted to simulate 
gaseous nebulae around ionizing stars, from the hot and ionized 
region to the neutral and cold gas. A grid of models was obtained 
for stellar and nebular parameters typical of planetary nebulae. 

The spatial distribution of molecules and the main mecha- 
nisms of molecular formation and destruction were discussed. 
The effect of the input parameters on the chemical composition 
was also analyzed. 

Our results show that a strong X-ray radiation field is needed 
to explain most of the column density ratios inferred from the ob- 
servations for species such as CO, CN, and HCO + . The molecu- 



lar chemistry is also sensitive to the dust-to-gas mass ratio, while 
the dependence on L* and is comparatively weak. 

Column density ratios involving HCN and HNC are, how- 
ever, not reproduced by the models. This problem is also found 
in other models of the literature. The HCN molecule is underes- 
timated in our models. We suggest that molecular data related to 
the reactions should be reviewed. An additional heating mecha- 
nism in the neutral region could be an alternative explanation. 

We showed that the A/(H2)/A/(CO) column density ratio de- 
pends on the value of N{CO). This is due to the distribution of 
both molecules inside the nebula. Using the relation obtained in 
the present work, we recalculated the molecular masses previ- 
ously obtained in the literature, showing that these masses are 
usually highly underestimated. As a result, we conclude that the 
missing mass problem can be solved by a more accurate estimate 
ofAT(H 2 )//V(CO). 
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